BCS-BEC crossover in spin-orbit coupled two-dimensional Fermi gases 
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The recent experimental realization of spin-orbit coupling for ultra-cold atoms has generated 
much interest in the physics of spin-orbit coupled degenerate Fermi gases. Although recently the 
BCS-BEC crossover in three-dimensional (3D) spin-orbit coupled Fermi gases has been intensively 
studied, the corresponding two-dimensional (2D) crossover physics has remained unexplored. In 
this paper, we investigate, both numerically and analytically, the BCS-BEC crossover physics in 2D 
degenerate Fermi gases in the presence of a Rashba type of spin-orbit coupling. We derive the mean 
field gap and atom number equations suitable for the 2D spin-orbit coupled Fermi gases and solve 
them numerically and self-consistently, from which the dependence of the ground state properties 
(chemical potential, superfluid pairing gap, ground state energy per atom) on the system parameters 
(e.g., binding energy, spin-orbit coupling strength) is obtained. Furthermore, we derive analytic 
expressions for these ground state quantities, which agree well with our numerical results within a 
broad parameter region. Such analytic expressions also agree qualitatively with previous numerical 
results for the 3D spin-orbit coupled Fermi gases, where analytic results are lacked. We show that 
with an increasing SOC strength, the chemical potential is shifted by a constant determined by the 
SOC strength. The superfluid pairing gap is enhanced significantly in the BCS limit for strong SOC, 
but only increases slightly in the BEC limit. 

PACS numbers: 03.75.Ss, 05.30. Fk, 74.20.Fg 



I. INTRODUCTION 

Spin-orbit coupling (SOC), the interaction between the 
spin and orbital degrees of freedom of a particle, has 
played important roles in condensed matter as well as 
atomic and nuclear physics. For instance, it is known 
that the coupling between electron spin and its linear 
momentum in solids leads to many important condensed 
matter phenomena such as spin and anomalous Hall ef- 
fects il|, i2|] , topological insulators and and superconduc- 
tors [3|, spintronics [1], etc. In atomic physics, the cou- 
pling between the electron spin and its motion around an 
atomic nucleus is responsible for many of the details of 
atomic structure @. 

Superfluidity and superconductivity are another im- 
portant phenomena in physics and have been widely 
studied in many physical systems, including solids. He- 
lium liquids, as well as ultra-cold atomic gases [6!]. Al- 
though particles in superfluids and superconductors usu- 
ally possess spins or pseudospins, the effects of SOC 
on superfluidity and superconductivity have remained 
largely unexplored. In this context, the recent experi- 
mental realization of SOC for ultra-cold atoms Q pro- 
vides a completely new platform for exploring many-body 
phenomena in spin-orbit coupled superfluids, including 
both Bose-Einstein condensates (BECs) j8j,|9|] and degen- 
erate Fermi gases [l3,[ll|- In the presence of SOC, vari- 
ous new and exotic superfluid phenomena may exist since 
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spins are not conserved during their motion. In particu- 
lar, in spin-orbit coupled BECs, the ground state phase 
diagrams as well as the collective excitations have been 
studied [3]. In spin-orbit coupled degenerate Fermi gases, 
the crossover physics from the Bardeen-Cooper-Schrieffer 
(BCS) superfluids of loosely bounded Cooper pairs to 
the BEC of tightly bounded molecules has been inves- 
tigated extensively, in both uniform and trapped three- 
dimensional (3D) gases [Ifl, [li|. However, the study of 
spin-orbit coupled two-dimensional (2D) Fermi gases is 
still lacked. 

On the other hand, the 2D degenerate Fermi gas (with- 
out SOC) in itself is one of the most active topics in 
ultra-cold atomic physics [l^]- Experimentally, the 2D 
degenerate Fermi cold atomic gases have been realized 
recently using a highly anisotropic pancake-shaped po- 
tential [13| , and the interaction energy in this system has 
been measured using the radio-frequency spectroscopy 
[l4| . In this system, tunable interactions between atoms 
through Feshbach resonance [l5| allow the exploration of 
the crossover physics from the BCS superfluids to the 
BEC of molecules in 2D ]]m\ . In addition to the study of 
many-body physics such as quantum fluctuations, non- 
Fermi liquid behavior, etc. [Tjl, the 2D Fermionic cold 
atomic gases are particularly interesting because of the 
existence of exotic topological excitations such as Majo- 
rana fermions with non-Abelian exchange statistics 18]. 

Motivated by these recent experimental breakthroughs 
in the realization of SOC and 2D Fermi gases, in this 
paper, we investigate, both numerically and analytically, 
the BCS-BEC crossover physics in 2D degenerate Fermi 
gases in the presence of a Rashba type of SOC. Our main 



results are the following: 

1) We derive the mean field gap and atom number 
equations suitable for the 2D spin-orbit coupled Fermi 
gases. 

2) These two equations are solved numerically and self- 
consistently, from which the dependence of the ground 
state properties (chemical potential, superfluid pairing 
gap, ground state energy per atom) on the system param- 
eters (e.g., binding energy, SOC strength) is obtained. 

3) We derive analytic expressions for these ground 
state physical quantities, which agree well with our nu- 
merical results within a broad parameter region. Such 
analytic expressions also agree qualitatively with previ- 
ous numerical results for the 3D spin-orbit coupled Fermi 
gases [ifll, where analytic results are lacked. 

4) We find that with an increasing SOC strength a, the 
chemical potential is shifted by a constant —ma^ deter- 
mined by a. The superfiuid pairing gap and the ground 
state energy per atom are affected at the order of a^ 
rather than a^, which means that weak SOC does not 
affect the pairing gap and the ground state energy per 
atom. In the strong SOC regime, the pairing gap and 
the ground state energy are enhanced significantly in the 
BCS limit, but only increase slightly in the BEC limit. 
Although these analytic results are obtained for the 2D 
Fermi gases, they also provide qualitative understanding 
for the numerical results in 3D Fermi gases [l3|, where 
similar changes of the chemical potential and superfluid 
pairing order with respect to the SOC strength are ob- 
served but analytic results are lacked. 

The paper is organized as follows. Section II describes 
the physical system: the spin-orbit coupled 2D degen- 
erate Fermi gases, and the corresponding Hamiltonian. 
In section III, we derive the mean field gap and atom 
number equations. These equations are self-consistently 
solved in section IV, both numerically and analytically 
(through perturbative methods), to obtain the ground 
state properties (chemical potential, superfiuid order pa- 
rameter, and ground state energy per atom) of the spin- 
orbit coupled Fermi gases in the BCS-BEC crossover. 
Section V consists of discussion and conclusion. 



II. THE PHYSICAL SYSTEM AND THE 
HAMILTONIAN 

The physical system in consideration is a 2D de- 
generate Fermi gas. Experimentally, the 2D degener- 
ate Fermi gas has been realized using a ID deep op- 
tical lattice along the third dimension, where the tun- 
neling between different layers is suppressed completely 
[mill. The ID optical lattice potential Vq sin^(27rz/Au,) 
can be generated using two counter-propagating laser 
beams (parallel to the z axis with a wavelength A^,). 
In this case, the two-body binding energy is given 
by E), ~ 0.915fiuJLGxp{y2TTl]^/as)/TT, where ujl = 
^8'!T^Vo/{mX'^) is the effective trapping frequency along 
the z axis, II = \/nJ(niujL), and Os is the 3D s-wave 



scattering length [10|. Therefore the two-body binding 
energy Ef, can be tuned by varying the s-wave scattering 
length Qs via the Feshbash resonance for the study of the 
BCS-BEC crossover physics flGj. For a small attractive 
fls — T' 0~, i?b — >■ 0, corresponding to the BCS limit. While 
for a small repulsive a^ — >■ 0^, E'b — >■ oo, corresponding 
to the BEC limit. When Ei, increases from to cx), the 
system evolves continuously from a BCS superfiuid to a 
BEC of molecules. 

The SOC for cold atoms can be generated by the in- 
teraction between atoms and laser beams, as shown in 
many previous literatures [20|, and demonstrated in a 
recent benchmark experiment [7|. In this paper, we con- 
sider only a Rashba type of SOC, and the effects of other 
types of SOC (e.g. Dresselhaus or the combination of 
both) can be investigated similarly. For simplicity we 
consider a uniform Fermi gas and neglect the weak har- 
monic trap in the 2D plane, whose effects can be incor- 
porated using the local density approximation. 

The Hamiltonian for this uniform 2D spin-orbit cou- 
pled degenerate Fermi gases can be written as 



where 



H — Hp + Hi + Hsoc, 



Hf - 2jCkCk<TCko 



(1) 

(2) 
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is the single atom Hamiltonian, C^^ is the creation oper- 
ator for a Fermi atom with the momentum k = {kx, ky), 
a =t, i are the pseudospins of atoms. Ck = ^k — A* with 
the kinetic energy ek = k'^ /2m, the chemical potential /i, 
and the atom mass m. Henceforth, we take the Planck 
constant h = 1. The s-wave scattering interaction be- 
tween atom 



Hi - ff X! '^-kt^ki^'kiC-kf, 



(3) 



where g is the effective scattering interaction parameter. 
In a 2D Fermi gas, 



^{2e^ + E,) 
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(4) 



The Hamiltonian for the Rashba type of SOC for atoms 
can be written as 

Hsoc = aY^[{ky + ikx)Cl^Cki + {ky-ikx)Cl^Cut], (5) 



where a is the SOC strength. In solid state materials, a 
is generally much smaller than Kp/ {2m) with Kp as the 
Fermi vector. However, in ultra-cold neutral atoms, a 
can reach the order of Kp/{2ra) ^2(/|. Such strong SOC, 
together with tunable interactions through the Feshbach 
resonance, may yield some exotic many-body phenom- 
ena that have not been explored in solid state systems. 
Recently, a generalized SOC with the Rashba and the 
Dresselhaus terms in the Hamiltonian ^ has been con- 
sidered Inl. 



III. MEAN FIELD GAP AND ATOM NUMBER 
EQUATIONS 

As the first step for the eventual understanding of the 
2D spin-orbit coupled Fermi gas, we consider the zero 
temperature superfluid physics under the mean field ap- 
proximation [l6|. Generally, the mean field approxima- 
tion call give qualitatively but not quantitatively correct 
results [2l| . In the mean field approximation, the super- 
fluid order parameter is taken as 



k 



(6) 



With this pairing order parameter, we can rewrite the 
two-body interaction Hamiltonian ^ as 

H, = -AV.9 + A 5](CHC_kt + Cl^t^l). (7) 



Therefore the total Hamiltonian can be rewritten as 

1 ^ _. . A2 — 
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H, 



B 



-^vI/t(k)Afk*(k) 



k 



Ck 



(8) 



under the Nambu spinor basis ^(k) = 

(Ck-f, C'kj., C'lki ; ^C'Lj^^)^, where Bogoliubov-de-Gennes 
operator 

Ck ak+ A 

afc_ Ck A 

A -Ck -ak+ 

A -ak_ -Ck 



Mk 



(9) 



preserves the particle-hole symmetry, k± = ky ± ik^- 
The quasiparticle excitation spectrum 



-B, 



k,± 



= AV(ek-M±afc)' + A2 



(10) 

is the eigenvalue of the matrix Afk, A = ± correspond 
to the particle and hole branches of the spectrum. For 
each branch, there are two different excitations due to 
the existence of the SOC. From Eq. ([5]), we see the total 
ground-state energy is 

A2 , 



E, 



1 



G 



^-^\l 



[C^^-7;{E+ + 



K- 



)]■ 



(11) 



Without the SOC, the term {E+_^_ 
the well-known form 



E^ _)/2 reduces to 



Eu 



v/(ek-M)' + A2 



(12) 



in the BCS theory. 

The ground-state properties of the 2D spin-orbit cou- 
pled Fermi gases can be obtained from the atom number 
equation 



dEc 



and the superfluid gap equation 



2^ d^i 






)], (13) 



dE, 



G 



l,dE, 



k+ 



dE, 



k- 



aA 



k " 



■ Eh 



OA 



dA 



)]=0. (14) 



IV. GROUND STATE PROPERTIES 
A. Numerical results 

We numerically solve the above atom number equa- 
tion ([T^ and the gap equation ([H)) self-consistently to 
obtain various ground state quantities. In Fig. 1, we 
plot the dependence of the ground state quantities: the 
chemical potential fj,, the superfluid order parameter A, 
and the ground state energy per atom E = Ec/n, on 
the physical parameters: the SOC strength a and the 
binding energy Eb- We see that the chemical potential /x 
decreases with the increasing spin-orbit coupling strength 
a. With increasing binding energy, the chemical poten- 
tial also decreases, signaling the crossover physics from 
the BCS superfluids to the BEC molecules. One inter- 
esting feature shown in Fig. I(a2) is that the shift of the 
chemical potential induced by the SOC depends only on 
the SOC strength. 

In Fig. l(bl), we see the superfluid order parameter 
A increases with increasing a. Here Aq is the superfluid 
order parameter without SOC. For a small a, the change 
of A/Aq is very small. The growth of A becomes signifi- 
cant only when aKp is large than the Fermi energy Ep. 
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Figure 1: (Color online) (al-cl) Plot of the chemical poten- 
tial fi, the dimensionless superfluid pairing gap A/Ao, and 
the ground-state energy per atom E with respect to the SOC 
strength a for the two-body binding energy Ei, — O.lEp 
(Black line) and 1.0 Ep (Red line). (a2-c2) Plot of fx, A/Ao, 
and E with respect to Ei, for aKp = 0.5Ep (Black line) and 
1.0 Ep (Red line). In all figures, the solid lines represent the 
approximate analytical results obtained in the paper, while 
the open symbols correspond to the exact numerical results. 



Furthermore, the effects of the SOC are clearly seen for 
a small Ef, (the BCS side), but not important for a large 
Eb (the BEC side). 

The ground state energy per atom E has a similar be- 
havior as A/Aq. It increases with a, but the growth is 
only important for a large a. Without SOC E = —^Ep, 
and the SOC induces a small correction. The change of 
E is significant in the BCS side (small Ei,), but negligible 
in the BEC side (large £'{,). We note that similar changes 
of the chemical potential and the superfiuid order param- 
eter have been observed numerically in a 3D spin-orbit 
coupled Fermi gases. All these numerical observed phe- 
nomena will be explained in the next subsection where 
the analytic expressions for these ground state quantities 
are derived. 



B. Analytic results 

Although the above numerical results give the depen- 
dence of the ground state quantities on the SOC strength 
for certain parameters, analytic results are desired for a 
better understanding of the underlying physics. In the 
following, we present a perturbative method (with a as 
the small parameter) to analytically extract the funda- 
mental ground-state properties of the spin-orbit coupled 
Fermi gases. For this purpose, we rearrange Eq. (jlip 
into two parts 



with 



Eg = Eq + Es, 



^°--^+E> 



Ek 



(15) 



(16) 



as the ground-state energy without SOC. For a 2D Fermi 
gas, Eq can be obtained exactly, yielding [16[ 



E. ^ £[A-ln( ^^-^^ )-My^^^TA^+.)-^]. 



A2 

2 

(17) 



The second term in Eq. ()15|) is given by 

i?soc=Ej^k-^(i?^,+ +i?^,J], (18) 

which describes the contribution from the SOC. 

Because it is difficult to derive a simple analytic ex- 
pression for ii^soc directly, we first perform a Taylor ex- 
pansion with respect to the SOC strength, and then do 
the summation over k. Formally, -Egoc can be written as 



first six orders 



A2- 



m), 



-4(vV 

3^ ^Vm' + A2" 
-16A2 

15(^^2 -1- A2)3/2 ' 

32/iA2 

35(a*2 + A2)5/2 ' 
64A2(A2 ~ V) 

315(Ai2 + A2)7/2' 

-128AiA2(3A2 - 4^2 
693(^^2 + A2)9/2 



(20) 
(21) 

(22) 
(23) 
(24) 
(25) 



Although the expression for Esoc is very complicate, 
the expressions for the chemical potential /i and the su- 
perfiuid pairing order A are very simple, as we will show 
later in the paper. With the ground-state energy Eq, 
the superfiuid pair order and the chemical potential can 
be derived by self-consistently solving the corresponding 
gap and number equations 



ln( 



^^i^ + A^-fi 
Eb ' 

v/^^ + A2 



M 



Ga(?7,A,^), (26) 

2Sp-G^(r,,A,Ai) (27) 



'-^W^V\ and 



analytically, where Ga(?7, A, /x) = ^ 

G^(?7,A, /x) = Yji ^^'l^'^'^'> r]\ Ep is the Fermi energy 

for a 2D non-interacting Fermi gas without SOC and 

with the density n = mEp/TT. Without SOC (77 — 0), 

Ga(??,A,^) = Gi,{tj,A,^) = 0, and Eqs. 1^ and ^ 

become 



Va^ 



A2 



vJ^ 



A2 



M = Eb, 
/i = 2Ep. 



In this case, the superfiuid pairing order and the chemical 
potential aregiven exactly by Ao = y/2EbEp and /io ~ 
Ep — Eb/2 [1^1, as expected. 

In the presence of SOC {rj 7^ 0), the nonlinear equa- 
tions (pS)) and (P7)) cannot be solved exactly. However, 
approximate solutions can be derived for the physical 
parameters within current experimentally achievable re- 
gion. Since the ground-state energy depends on A2, the 
solutions of Eqs. (I^5|) and (l?7l) can be assumed to be 
fi = ^i^of^iV^ and A2 = X)i=o^»^*- Substituting these 
expressions into the nonlinear equations ([^5]) and ^7} 
and then comparing the coefficients for the same order of 
77*, we obtain A^ and /i^ respectively. 



E. 



SOC - ■^C^{^l, A)!]"- 



(19) 



with 77 — mo? 12. The coefficients Gi(/x, A) can be ob- 
tained, in principle, for any order. Here, we only give the 



1. Chemical •potential 

After a straightforward but tedious calculation, the co- 
efficients for the chemical potential are given by pb\ = 



-2, 1^2 = ^Eb/[3{Eb + 2Ef)^] and ^3 = -64EbiEb - 
4:Ep)/[15{Eb + 2Ef)^]- Note that the second order ^2 
is already very small for all different values of Ef, when 
77 < Ep, therefore the high-order terms do not affect the 
chemical potential significantly. The chemical potential 
can then be written as 



Hc,Ep~- 



2?7. 



(28) 



From Fig. la, we see Eq. ((28|) agrees well with the exact 
values of the chemical potential obtained from numeri- 
cally solving Eqs. P^ and (|14p self-consistently. 

Equation ((28|) shows clearly that with the increasing 
SOC strength a, the chemical potential /i is decreased 
by 2?7 = TTia^, as shown in Fig. 2. If we define an effec- 
tive chemical potential 'p = fi + 2r] = fi + ma^ , Eq. (pS)) 
can be rewritten asjl — Ep — £'b/2. Similar to the previ- 
ous discussion without SOC, wc find that in the asymp- 
totic BCS limit with a weak bound state [Ei, <^ Ep), 
the chemical potential JI ~ Ep. However, in the deep 
BEC regime with a strong bound state {Ef, ^ Ep), the 
chemical potential Jl = —Eb/2. The BCS-BEC crossover 
may occur around the crossover point /I = [16| . Fur- 
thermore, in the presence of SOC, Eq. ([^5)) can also 
be written as i-i — fif + /ia, where /ij = Ep — 2ry 
is the chemical potential for the free Fermi gas and 
/lA = —Eb/2 + ^i^i fiAi{Eb)T]^ reflects the revision of 
the chemical potential induced by the two-body binding 
energy Eb- Here /^ai = implies that the binding energy 
Eb has no effect on the chemical potential /x at the order 



of rj 



72. 



It is important to point out that although the chemi- 
cal potential fi only has a simple shift from that without 
SOC, the underlying physics is quite different. Without 
SOC, the pairing wave function is simply singlet. How- 
ever, in the presence of SOC, each energy band contains 
both spin up and down components. As a result, the 
pairing wave function has a more complicated structure 
with both singlet and triplet components 22^ . The triplet 
pairing correlations in s-wave superfluids may be used 
to detect the anisotropic Fulde-Ferrell-Larkin-Ovchinikov 
states [231 • In addition, the 2D Fermi gases with SOC can 




Figure 2: (Color online) (a) The chemical potential without 
SOC in the BCS limit, (b) The chemical potential with SOC 
in the BCS limit. For weak SOC, their difference A^ = ma^. 



exhibit p-wave character in the helicity bases. If a Zee- 
man field in Hamiltonian ([8]) is added, a novel topolog- 
ical phase transition from non-topological superfluids to 
topological superfluids can be induced. In the topological 
superfluid phase, there exist Majorana fermions and the 
associated non-Abelian statistics, which are the critical 
ingredients for implementing topological quantum com- 
puting [18,] . 



2. Superfluid order gaps 

The effects of SOC are more interesting for the super- 
fluid pairing gap A. Applying the same procedure as 
that for the chemical potential, we find 



A^ ~ 2EbEp 



WEbEp 
i{Eb + 2Ep) 



(29) 



There is no first-order correction with respect to -q {^ a^) 
for A^, and the second-order rj^ (~ a^) is the leading cor- 
rection. Moreover, the second-order coefficient d"^ t^ jdrf' 
is always positive and has a maximum 4?7^/3 (= rn^a^/S) 
when Eb = 2Ep. The high-order coefficients for rf and 
774 arc A3 = -h\2EbEp{Eb - Ep)l\Vh{Eb + 2Ep)^\ and 
A4 = %\EbEpl([S)YJEl-Zh^%EbEp^\^l%E'],)l\i\^{Eb^ 
2Ep)^\ respectively. 

In order to see the effect of SOC on the superffuid 
pairing gap more clearly, we introduce a dimensionless 
quantity 



Arf = 




7>{Eb + 2EpY 



(30) 



where Aq is the superfluid pairing gap without SOC. 
In the asymptotic BCS limit with a weak bound state 
{Eb < Ep), Ad ~ yrrVTeJ. For a weak SOC 
{rj <C Ep) in typical solid-state materials, A^ ~ 1, 
which means that the SOC does not affect the super- 
fluid pairing gap significantly. However, for a strong 
SOC r] ^ Ep that has been achieved for ultracold atoms 
[3], this gap can be enhanced greatly. In the deep BEC 
regime with a strong binding energy (Eb 3> Ep and 
Eb :^ T]), Ad ~ yTT8rfJ{3ElJ , therefore the super- 
fluid pairing gap increases only slightly with the increas- 
ing SOC strength. These analytic results agree well with 
the numerical results shown in Fig. lb. Note that similar 
behavior for A^ is also observed in the numerical results 
in 3D M. 



3. Ground-state energy per atom 



In terms of Eqs. p8p and ([29| . the ground-state energy 
per Fermi atom E = Eq/ti can be obtained 



E: 



1 

-Ep 



8Ef 



HEb 



2Ep)^'^ 



(31) 



The comparison of Eq. (|3T|) with the direct numerical 
simulation results is shown in Fig. Ic. The ground- 
state energy, like the superfluid pairing gap, depends on 
rf (^ a^). In the asymptotic BCS Hmit {Ei, <^ Ep), 
E ~ ~i?i?/2 + 277^/(3i?F), which means that the ground- 
state energy can be enhanced significantly only for a 
large 77. In the deep BEC regime {Ei, ^ Ep and 
Eb > 7]), E ~ -Ef/2 + SEprj^/iSEl), therefore the 
ground-state energy only increases slightly even when 
77 ~ Ep. In addition, the high-order coefficients for 77"^ 
and rj'^ are given by E3 = -l28EbEp/[5{Ei, + 2Epy] 
and E4 = 128EbEp{A91Eb - 538 E p) / [315 {Eb + 2Epf] 
respectively. Note that although such mean field the- 
ory gives qualitatively correct results, it may not agree 
quantitatively with the experimental results in the deep 
BEC regime as shown in recent Monte Carlo numerical 
simulation for the 2D Fermi gases without SOC |2l| . 

Finally, we want to remark that if the high-order terms 
in the chemical potential /i, the superfluid pairing gap A, 
and the ground-state energy per Fermi atom E — Ea/n 
are included, the analytical results in Eqs. (pS)) . ([50)) 
and (j3ip fit better with the numerical results even for 
Ep < aKp < 3.0Ef. 



V. DISCUSSION AND CONCLUSION 

The mean field zero temperature BCS-BEC crossover 
physics discussed above provides the first critical step 
towards the understanding of the 2D spin-orbit coupled 
degenerate Fermi gases. Clearly, many issues need be 
further explored in the future, as demonstrated by the 
development along this direction after the initial submis- 
sion of our paper. In particular, the finite temperature ef- 
fect need be taken into account in a realistic experiment. 
Without SOC, it is known that there are no superflu- 
ids for 2D degenerate Fermi gas at a finite temperature 



24|, where relevant physics is the Berezinskii-Kosterlitz- 
Thouless (BKT) transition [23 , leading to the generation 
of the vortex-antivortex pairs. Recent some interesting 
effects of the SOC on the BKT transition 26] has been 
investigated. In experiments, the Fermi gases are con- 
fined in a 2D harmonic trap, whose effects may be taken 
into account using the local density approximation, as 
shown in the recent preprint [26|. The trap geometry 
for the Fermi gases may also be quais-2D instead of the 
strict 2D. In this case, the confinement along the third 
direction may affect the critical transition temperature, 
which need be further explored. 

In summary, motivated by the recent experimental 
breakthrough for the realization of the SOC for cold 
atoms and the 2D degenerate Fermi gas, we investigate 
the zero temperature BCS-BEC crossover physics in 2D 
spin-orbit coupled degenerate Fermi gases using the mean 
field approximation. By solving the corresponding gap 
and atom number equations both numerically and an- 
alytically, we reveal the ground state properties of the 
spin-orbit coupled 2D Fermi gases. Our analytic results 
agree quantitatively with our numerical results in 2D and 
qualitatively with previous numerical results for the 3D 
spin-orbit coupled Fermi gases, where analytic results 
are lacked. The analytic expressions for various physi- 
cal quantities may provide a powerful tool for engineer- 
ing and probing many new topological phenomena in 2D 
Fermi gases, including the intriguing Majorana physics 
and the associated non-Abelian statistics. 
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